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Abstract:  Fermi  y  -ray  telescope  data  of  GRB  080916C  with  ~  1055  erg  in  apparent  isotropic  y  -ray  energy,  show  a 
several  second  delay  between  the  rise  of  100  MeV  -  GeV  radiation  compared  with  keV  -  MeV  radiation.  Here  we  show 
that  synchrotron  radiation  from  cosmic  ray  protons  accelerated  in  GRBs,  delayed  by  the  proton  synchrotron  cooling 
timescale  in  a  jet  of  magnetically -dominated  shocked  plasma  moving  at  highly  relativistic  speeds  with  bulk  Lorentz  factor 
r  ~  500 ,  could  explain  this  result.  A  second  generation  electron  synchrotron  component  from  attenuated  proton 
synchrotron  radiation  makes  enhanced  soft  X-ray  to  MeV  y  -ray  emission.  Long  GRBs  with  narrow,  energetic  jets 

accelerating  particles  to  ultra-high  energies  could  explain  the  Auger  observations  of  UHE  cosmic  rays  from  sources 
within  100  Mpc  for  nano-Gauss  intergalactic  magnetic  fields.  The  total  energy  requirements  in  a  proton  synchrotron 
model  are  °c  rI6/3 .  This  model  for  GRB  080916C  is  only  plausible  if  T<500  and  the  jet  opening  angle  is  ~1  . 
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1.  INTRODUCTION 

An  integrated  fluence  of  2.4  xlCT4  erg  cm-2  was 
measured  from  GRB  080916C  with  the  Large  Area 
Telescope  (LAT)  and  Gamma  ray  Burst  Monitor  (GBM)  on 
the  Fermi  Gamma  ray  Space  Telescope,  with  one  third  of 
this  energy  in  the  LAT  [1],  At  a  redshift  z  =  4.35 ±0.15  [2], 
GRB  0809 16C  has  the  largest  apparent  energy  release  yet 
observed  from  a  GRB.  A  significant  ~  4.5  s  delay  between 
the  onset  of  >100  MeV  compared  to  the  ~  8  keV  -  5  MeV 
radiation  is  found  (the  characteristic  duration  of  the  GBM 
emission  is  »  50  s). 

The  spectrum  of  GRB  080916C  was  fit  by  the  smoothly 
connected  double  power-law  Band  function  [3]  to  multi-GeV 
energies,  though  with  changing  Band  spectral  parameters 
and  peak  photon  energy  in  different  time  intervals.  The 
emergence  of  delayed  spectral  hardening  is  represented  by  a 
Band  beta  spectral  index  changing  from  f>  =  -2.6  in  the  first 
3.6  seconds  following  the  GRB  trigger  to  /3  =  -2.2  at  later 
times  [1],  Here  we  show  that  a  hard  spectral  component 
arising  from  cosmic-ray  proton  synchrotron  radiation 
explains  the  delayed  onset  of  the  LAT  emission.  If  GRBs 
accelerate  UHECRs,  then  the  delayed  onset  of  the  LAT 
emission  after  the  GBM  trigger  should  be  a  regular  feature  of 
GRB  spectral  evolution. 

2.  PROTON  ACCELERATION  AND  RADIATION 

GRB  blast  wave  calculations  usually  treat  electrons  [4, 
5],  but  protons  and  ions  will  also  be  accelerated  if  they  are 
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present  in  the  relativistic  flows  in  black-hole  jet  systems. 
Here  we  consider  protons  accelerated  in  GRB  blast  waves  to 
such  energies  that  they  can  efficiently  radiate  hard  ~  GeV  - 
TeV  photons  by  the  proton  synchrotron  mechanism  [6-11]. 
The  highest  energy  photons  are  reprocessed  by  yy  —*  e+e~ 
opacity  [12]  to  make  an  injection  source  of  electrons  and 
positrons  that  cool  by  emitting  <  GeV  electron  synchrotron 
radiation.  Two  delays  arise,  the  first  from  the  time  it  takes  to 
accelerate  protons  to  a  saturation  Lorentz  factor  where  the 
acceleration  rate  equals  the  synchrotron  rate.  A  second  delay 
arises  during  which  sufficient  time  passes  to  build  up  the 
spectrum  of  the  primary  protons  so  that  they  become 
radiatively  efficient  in  the  LAT  band.  Observations  of  GRB 
0809 16C  are  consistent  with  the  delayed  onset  between 
GBM  and  LAT  emission  being  caused  by  the  second  delay 
where  the  evolving  proton  cooling  synchrotron  spectrum 
sweeps  from  higher  energies  into  the  LAT  waveband.  This 
radiation  is  emitted  from  a  jet  of  magnetically-dominated 
shocked  plasma  with  'QK  >  I  ,  where  tH  is  the  ratio  of 
magnetic-field  to  proton/particle  energy  density  in  the 
plasma. 

The  rapid  variability,  large  apparent  luminosity,  and 
detection  of  high-energy  photons  from  GRBs  can 
be  understood  if  this  radiation  is  emitted  from  jetted 
plasma  moving  with  bulk  Lorentz  factor  T»1  towards  us. 
Detection  of  3  GeV  and  13  GeV  photons  from  GRB 
0809 16C  suggests  T3=r/1000~l  (Ref.  [1]  and  below, 
Section  3).  For  variability  times  tv  ~  I  s  and  instantaneous 
energy  fluxes  <1>  =  10~3<!l5  erg  cm-2  s’1,  the  internal 
radiation  energy  density  in  the  fluid  is 
ur  =  4jrdl$>  /  (4jtR2lT2  )  «  (1  +  z?  d2L<S>  /  (TV/2 ) ,  implying 
a  characteristic  jet  magnetic  field  of 
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B  (kG)  ~2^°Pf-5  -  2gB2  pb2t~1'2  (s),  E~m  (1 3GeV)  (1) 

r3fv(s) 

where  primes  refer  to  the  comoving  frame,  pb  is  the  baryon- 
loading  parameter  giving  the  relative  energy  in  nonthermal 
protons  compared  to  y  -rays.  The  last  relation  in  equation 

(1)  assumes  T « rmjn  from  the  opacity  condition  t  =1, 
which  can  be  written  as 


r  s 

min 


oTdl(l  +  z)2  fisl 


4  tjn/y 


(2) 


where  f-  is  the  vFv  flux  at  photon  energy 
e  =  2T2 /[(l  +  z)2e1]  (in  mec 2  units)  and  E  =  mec2£l  is  the 
highest  energy  photon.  Thus  the  total  jet  energy  is 
°c  gBphl I>  .  A  characteristic  field  B  ~  10  -  100  kG  is 


consistent  with  the  absence  of  a  distinct  self-Compton 
synchrotron  component  in  GRB  0809 16C  [1],  which 
happens  when  the  magnetic  energy  density  is  larger  than  the 
nonthermal  electron  energy  density,  or  tBPb  >>  1  ■ 


In  Fermi  acceleration  scenarios,  protons  gain  energy  on 
timescales  exceeding  the  Larmor  timescale,  implying  an 
acceleration  rate  yaccp  =  (p~'eB  /  mpc  ,  where  0“ 1  <<  1  is  the 

acceleration  efficiency.  Equating  the  acceleration  rate  with 
the  synchrotron  loss  rate  gives  a  saturation  Lorentz  factor  for 
protons,  namely 


■  2xi°8 
sa,’P  me{<pB’ )  4af  ^/(0/lO)fi5  ’ 
where  =1/137  is  the  fine  structure 


(3) 

constant. 


Bcr  =  m2c 3  /  eh  =  4.4  x  1013  G  is  the  critical  magnetic  field, 
and  the  mean  magnetic  field  of  the  radiating  region  is 
B  =  I  05  B5  G.  The  observer  measures  a  time 


=  l+_ z'n}c_  I  6 7i<j>  0.01^/0/10 

I  me  Y  earB  3  ~  I  jC  '  S 
for  protons  to  reach  ysat  . 

The  proton  synchrotron  saturation  frequency  (in  mec 2 
units),  corresponding  to  the  proton  synchrotron  frequency  of 
protons  with  yp  =  ysat  p ,  is 


1 +z  m.  8 a. 


(5) 


m'p  1  +  z  mup  1  +  z  me  8  af  0/10 

This  frequency  corresponds  to  a  photon  energy  of  »  8 
TeV  for  GRB  0809 16C.  Analogous  to  electron  blast  wave 
physics  [13],  proton  synchrotron  losses  make  a  cooling  break 
at  the  proton  cooling  Lorentz  factor  yc  =  ym  (fm  / 1) , 


which  is  obtained  by  equating  the  synchrotron  energy-loss 
timescale  with  the  comoving  time  t  =17 /(1  +  z). 
Consequently  the  proton  synchrotron  cooling  frequency  for 
protons  with  yp  « yc  p  is  ecp(t)  =  £ml  p(tsal/t)2  for  our 

simplified  model  of  continuous  acceleration  and  uniform 
injection  of  particles  during  the  first  ~  8  s  to  the  observer. 


Photons  with  energies  above  ~  1-10  GeV  are  strongly 
attenuated  through  yy  — *  e+e~  processes  in  the  source, 

inducing  a  nonthermal  e2  injection  that  makes  a  second- 
generation  electron  synchrotron  component.  Proton 
synchrotron  photons  radiated  by  protons  with  yp  =  ysal 
materialize  into  electrons  and  positrons  with  Lorentz  factor 
ye  «  £sat  /  2  .  The  electron  synchrotron  saturation  frequency 

from  this  second  generation  of  synchrotron  emission  is 
observed  at 


3  T  2  B  ( mp  27  V 
2 1  +  z  Bcrmel6afi 


io3r  3b5 

(0/1O)2  ‘ 


(6) 


Therefore  the  second-generation  synchrotron  radiation  cuts 
off  above  » 600  MeV  for  GRB  0809 16C.  This  will  be 
below  the  LAT  sensitivity  because  only  the  proton 
synchrotron  radiation  from  protons  with  yp  =  ysat 
contributes  to  the  second-generation  electrons  and  positrons 
that  make  emission  at  £  =  £sat  e  and  the  number  of  such 
protons  is  small  for  a  steep  injection  proton  spectrum. 

The  direct  proton  synchrotron  radiation  from  cooling 
protons  equals  the  electron  synchrotron  saturation  frequency 
when  £cp(t)  =  £mlp(tsa,  It)2  =  £sale ,  which  is  observed  to 
take  place  at  time 


LK  mt  64 a f  _41  +  z  mpcBcr 
v  B  mp  81  3  T  eB2 


(7) 


Lor  GRB  080916C,  td « 1.4(0/ 10)/ s,  which 
corresponds  (depending  precisely  on  £sal  e )  to  the  time 

required  for  proton  synchrotron  radiation  to  become  strong 
in  the  LAT  energy  band.  We  propose  this  effect  as  the  reason 
for  the  delayed  onset  of  the  LAT  emission. 

In  addition  to  the  direct  proton  synchrotron  radiation,  pair 
synchrotron  radiation  is  formed  from  the  internal  attenuation 
of  proton  synchrotron  photons  to  make  ultra-relativistic 
pairs.  If  protons  are  accelerated  with  a  number  index  k  (i.e., 
N(y)  oc  y~k ),  then  the  cooling  spectrum  breaks  from  a  proton 
distribution  with  index  k  to  a  steeper  one  with  index  k  + 1 . 
The  proton  synchrotron  vFv  flux,  in  the  absence  of  yy 
opacity,  is  simply 
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where  =  'Vp&GBM  -  y’min.p  -  r«-,  -  the  relative  Lorentz 

factor  of  the  relativistic  wind  in  the  jet  and  the  shell  material, 
and 


k-  2\\-{eminJemUprkV\ 

2  /  (£mm,P  /  £sol,p)C2_i)/2  “  1  ' 


(9) 


^  =  [ln<  Ym.p  1  Ynan.p  )]"‘  When  k  =  2. 

The  direct  proton  synchrotron  flux  from  GRB  0809 16C 
is  attenuated  at  all  energies  from  e  ~  I  -  10  GeV  to 

£  « £ sM  ,  equation  (5),  to  make  a  second-generation 

electron  injection  spectrum  with  the  same  form  as  equation 
(8),  though  with  subscript  p  —*  e  spectral  indices 
(3 - &) / 2  — »  (3 - £) / 4  and  (2-k)/2  -*(2-£)/4  and 

re.syn  _  j_  y-p.^n  seconcj-generation  spectrum  has  a  low- 

^  bsat,p 

energy  cut-off  related  to  £  ,  below  which  it  receives  no 

further  injection  pairs.  Note  that  £ce(f )  =  £sate(tsat  /  04 ,  where 
£lii(  e  is  given  by  equation  (6). 


Fig.  (1).  Synchrotron  model  for  the  y  -ray  spectrum  of  GRB 
080916C.  The  vFv  flux  in  the  GBM  range  is  produced  by  primary 

nonthermal  electron  synchrotron  radiation.  The  heavy  dark  curves 
show  the  fitted  Band  spectrum  for  intervals  (a),  from  0  to  3.6  s  after 
the  trigger,  and  interval  (b),  from  3.6  s  to  7.7  s.  A  strong  proton 
synchrotron  component  is  formed  at  high  energies  and  sweeps  into 
the  LAT  band  after  several  seconds,  thus  making  the  time  delay 
between  the  GBM  and  LAT  emissions.  A  weak  pair  synchrotron  y  - 
ray  component  is  fomied  in  the  GBM  and  LAT  band  from  second- 
generation  pair  synchrotron  radiation  formed  by  proton  synchrotron 
photons  that  are  attenuated  by  yy  processes  in  the  jet.  Here  we 
assume  uniform  injection  over  8  s,  and  take  <f>  =  10  .  At  a  fixed 
photon  energy  below  the  cooling  frequency,  the  vFv  flux  increases 
a  t  for  both  the  proton  synchrotron  and  pair  synchrotron  radiation. 


Fig.  (1)  shows  the  evolving  proton  synchrotron  spectrum 
for  constant  injection  with  time  using  the  above  relations  for 
parameters  of  GRB  0809 16C,  with  the  Band  function  in 
interval  (a),  the  first  3.6  s  following  the  trigger.  The  Band 
function  in  the  ~  10  keV  to  <100  MeV  range  can  be  due  to 
synchrotron  radiation  by  primary  electrons,  or  thermal 
radiation  from  a  jet  photosphere.  For  the  proton  spectrum, 
we  use  k  =  2.3  ,  pb=\l  ,  t,B  =  5  ,  and  ymin  -  10  .  In  this 

simple  picture  of  the  GRB  jet,  the  protons  are  accelerated  to 
the  saturation  Lorentz  factor  within  « 0.01  s,  making  a 
prompt  second-generation  electron  synchrotron  spectrum  too 
weak  to  be  detected,  followed  after  a  few  seconds  by  strong 
direct  proton  synchrotron  emission  that  sweeps  into  the  LAT 
band  from  high  energies.  The  time  delay  in  the  high-energy 
y- ray  flux  occurs  even  for  variable  injection,  provided  it 

occurs  over  ~  8  s  to  the  observer.  Although  the  data  for 
GRB  080916C  are  consistent  with  a  Band  function,  spectral 
analyses  with  the  addition  of  a  proton  synchrotron 
component  will  determine  whether  this  more  complicated 
model  is  compatible  with  Fermi  LAT  GRB  data. 

Implicit  in  the  treatment  is  that  the  wind  impacts  either  a 
shell  ejected  earlier  by  the  central  engine  or  some  roughly 
uniform  density  material  that  already  existed  in  the 
surrounding  medium.  The  extent  of  this  material  is 
(r)  =  cT2At /(1  +  z)  «  5  x  lO^T^Af  /8  s)  cm.  The  magnetic 
field  is  thus  assumed  to  be  roughly  constant  over  the  first 
~  10  s  because  of  the  constant  density  medium  that  is  being 
swept  up  at  the  shock.  This  assumption  is  supported  by  the 
roughly  constant  value  of  electron  synchrotron  energy  flux 
<1>G,(W  measured  with  the  GBM,  in  the  same  time  interval 
[1]. 

3.  TOTAL  ENERGY  OF  GRB  080916C 


The  comoving  synchrotron  cooling  timescale  of  an  ion 
with  atomic  mass  A  and  charge  Z  is  given  by 
tsyn  =(A'  /  Z4)(mp  /  me)3(6jrmec  /  aTB2y).  The  comoving 
peak  synchrotron  photon  energy  is 

B  =  (Z/  A)(me  / mp)By 2  /  Bcr .  Equating  the  observer  time 


:(l  +  z)f  /T  for  the  emission  to  be  radiated  at 


measured  energy  mec  e  =  fe'  /  (1  +  z)  =  1 00£1( 


implies 


comoving 


magnetic 


MeV 

field 


B  (G)  =  2.0  x  10  A  Z~  E~qq  t„„  (s)  and  an  isotropic  jet 
power,  dominated  by  magnetic-field  energy,  given  by 

R2cT2B' 


Lb  -  ■ 


=  2x10 


58 


rf^410/3L2'3(5) 


7-14/3^2/3 
A  ^100 


erg  s  (10) 


letting  the  blast-wave  radius  R  =  r2ct  /  ( 1  +  z)  ■  For  Fe 
( A  =  56,  Z  =  26  ),  the  power  requirements  are  reduced  by  a 
factor  «0.17.  Here,  however,  we  consider  only  proton 
synchrotron  radiation. 

Eq.  (10)  shows  that  LB  oc  r16/3  [11].  The  absolute  jet 
power  varies  as  the  square  of  jet  opening  angle  6 .  .  In  Ref. 
[1],  only  the  uncertainty  in  the  redshift  was  used  to  provide 
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uncertainty  in  Tmjn  using  simple  yy  opacity  arguments, 
Furthermore,  validity  of  the  cospatial  assumption  that  the 
soft  photons  are  found  in  the  same  region  as  the  hard 
photons  was  assumed  in  the  calculation  of  Tmjn .  We  use  a 
likelihood  ratio  test  to  calculate  the  chance  probability  to 
detect  a  photon  with  energy  greater  than  the  maximum 
measured  photon  energy  in  a  given  time  bin.  For  a  -2.2 
photon  number  spectrum,  this  test  gives  values  of 
r  =  (0.90, 0.82, 0.75 )Tmin  for  exponential  escape  and 
r  =  (0.79, 0.58, 0.36)r„,„  for  slab  or  spherical  escape  at  the 
(1  cr  ,  2  cr  ,  3  cr  )  levels,  respectively. 

If  Type  Ib,c  supernovae  are  the  progenitors  of  long 
GRBs  like  GRB  080916C,  then  0.  >  0.8°  [14],  Taking 
a  conservative  lower  limit  T3  =  0.5  gives  the  absolute 
energy  requirements  for  GRB  0809 16C  of 

£  =  3xlO53(r3/O.5)16/3(0. /ldeg)2(^„ /8s)5/3  erg  after 
integrating  eq.  (10)  over  time  and  multiplying  by  a  two-sided 
jet  beaming  factor  fb  as  1.5  x  lO-4^  / 1  deg)2 .  The  rotational 
energy  available  in  a  core  collapse  supernova  could  be  as 
large  as  «  5  x  1054  erg  for  a  10  Msol  core  [15]. 

Although  it  appears  to  be  a  coincidence  that  T  «  Tmin  in 
our  model,  which  is  required  because  of  excessive  powers 
when  r>  r„„„,  more  complicated  geometries  might  relax  the 
bulk  Lorentz  factor  requirement  further.  If  the  inner  engine 
makes  the  prompt  MeV  radiation  and  residual  shell 
collisions  at  larger  radii  make  LAT  y  -ray  photons,  then  T 
could  be  as  low  as  ~  300  [16].  In  this  case,  the  absolute 
energy  release  could  be  as  low  as  »  2  x  1O52(0;.  / 1  deg)2  erg. 

Even  though  the  model  in  Ref.  [16]  provides  a  separate 
explanation  for  the  delayed  onset  (and  predicts  that  the 
variability  timescale  of  the  >  100  MeV  radiation  is  longer 
than  the  keV/MeV  radiation),  lack  of  cospatiality  does  not 
guarantee  delayed  onset.  For  example,  photospheric 
emission  with  leptonic  emission  from  internal  shells  would 
neither  be  cospatial  nor  necessarily  exhibit  a  delayed  onset. 

The  jet  break  time  with  apparent  isotropic  energy  release 
«  2x  1057  ergs  is  tbr  =  0.3(9  j  / 1  deg)'c/3n0,/3  d.  The  jet  break 
would  have  taken  place  before  Swift  slewed,  at  =  7], +  17.0 
hr,  to  GRB  080916C  [17],  with  a  hard  electron  spectrum  to 
explain  the  shallow  X-ray  decay  observed  >105  s  after  the 
GRB  [18]. 

A  *3  s  and  »2  s  delayed  onset  of  >100  MeV  emission 
detected  in  another  long  GRB  090902B  [19]  and  in  the  short 
GRB  081024B  [20],  respectively,  could  be  explained  with 
proton  synchrotron  radiation  as  we  discussed  here,  although 
lack  of  redshift  measurement  for  GRB  081024B  prohibits 
calculation  of  the  bulk  Lorentz  factor  and  the  total  y- ray 

energy.  It  is,  however,  interesting  to  note  that  Fermi-LAT- 
detected  short  GRBs  emit  higher  fluence  in  the  LAT  range 
than  in  the  GBM  range  [20,  21]  which  could  be  due  to 
hadronic  emission  processes. 


4.  ULTRAHIGH-ENERGY  COSMIC  RAYS  FROM 
GRBS 

GRBs  have  long  been  considered  as  candidate  sources  to 
accelerate  UHECRs  [22,  23].  The  energy  of  protons  with 
y  =  ysal  ,  if  they  were  to  escape  from  the  GRB  blast  wave, 

is  =2xl02Or3  /f(j>n0)B5  eV,  so  GRB  0809 16C  can  in 
principle  accelerate  UHECRs.  The  Larmor  timescale 
tesc(s)  “  0.2  /  sj((j>  / 10)0/  at  y'p  =  y'sat  p  is  much  less  than  the 

light-crossing  time  l,c(s)  <=  200fv(s) ,  so  that  escape  depends 
on  transport  in  the  jet  plasma  magnetic  field  impedes 
escape  [24].  Photohadronic  processes  can  assist  escape 
by  converting  protons  to  neutrons,  but  are  more  important 
when  the  internal  photon  density  is  high,  or  when 
the  variability  timescale  tv  is  small  [25]. 

When  fv( s)  <  0.01(d>GBM  / 10“5  erg  cnrV1  )/(T43£pk) ,  then 

photopion  losses  are  important  in  GRB  0809 16C  (see  Ref. 
[26]  for  applications  to  GRB  090510).  The  shortest 
variability  timescale  for  GRB  0809 16C  observed  with 
INTEGRAL  is  0.1  s  [2],  so  unless  the  corresponding  size 
scale  for  the  radiating  region  was  even  shorter,  we  can 
neglect  photohadronic  processes.  But  a  shorter  variability 
timescale  would  require  a  larger  bulk  Lorentz  factor  T  in 
order  to  explain  detection  of  multi-GeV  photons,  in  which 
case  photohadronic  efficiency  is  reduced  [10]. 

We  estimate  the  rate  of  long-duration  GRBs  as  energetic 
as  GRB  080916C  within  the  =  100  Mpc  clustering  radius  for 
UHECRs  observed  with  Auger  [27,  28].  For  maximum  total 
energy  releases  of  «  1054  erg,  the  GRB  080916C  jet  opening 
angle  6j  <  100 /T  =  0.1/T3 .  The  inferred  GRB  rate  of 

=  2fh  Gpc  “3  yr  _1  [29]  at  the  typical  redshift  z  =  1  -  2  is  a 
factor  =1/10  smaller  at  =  100z/1GO  Mpc  due  to  the  star 
formation  rate  factor  and  a  factor  fb  >  200  larger  due  to  a 
beaming  factor.  A  60£60  EeV  UHECR  is  deflected  by  an 
angle  =  4.4°ZBnG£,^01^i11(,)oA11,2  in  intergalactic  magnetic  field 
with  mean  strength  BnG  nG  and  coherence  length  of  A, 
Mpc.  Deflection  causes  dispersion  in  time  of  arrival  of 
UHECRs  [30]  and  increases  the  apparent  rate.  The 
corresponding  number  of  GRB  sources  within  =  100  Mpc 
with  jets  pointing  within  4“  of  our  line-of-sight  is 
»  30 (fb  /  2OO)02GE(jo  A3/2  .  Complications  arising  from  non- 
uniform  magnetic  field  geometries  (e.g.  Refs.  [31,  32])  can 
lead  to  different  values  in  the  rate  estimates,  but  still  allow 
GRBs  as  UHECR  hosts.  Thus  if  typical  long  duration  GRBs 
have  a  narrow,  highly  relativistic  core  accelerating  UHECRs, 
then  long-duration  GRBs  could  account  for  the  Auger  events 
within  the  GZK  radius. 

5.  DISCUSSION 

We  have  developed  a  hadronic  model  based  on 
synchrotron  radiation  by  protons  in  a  highly  magnetized 
shocked  plasma,  to  explain  delayed  onset  of  high  energy  (> 
100  MeV)  emission  observed  with  Fermi  from  GRB 
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0809 16C.  This  proton  synchrotron  spectral  component,  in 
addition  to  the  electron  synchrotron  and  photospheric 
emission  [33]  that  dominates  in  the  keV  -  MeV  range, 
initially  starts  at  much  higher  energy  and  later  sweeps  into 
the  Fermi  LAT  range,  thus  causing  a  time  delay  in  the 
prompt  phase.  Our  model  is  compatible  with  an  internal 
shocks  scenario  and  could  also  be  consistent  with  an  origin 
of  the  high-energy  emission  in  GRB  0809 16C  from  an 
external  shock  [34].  The  late  time  extended  emission 
observed  in  GRB  080916C  and  other  GRBs  could  be  due  to 
the  slower  cooling  of  the  protons  in  the  early  afterglow 
related  to  the  external  shock  [6,  9]. 

Our  approach  differs  from  the  conclusions  of  [35]  that 
the  jet  outflow  energy  in  GRB  0809 16C  is  dominated  by 
Poynting  flux  rather  than  particle  energy.  Their  conclusion 
depends  on  the  assumption  that  the  engine  radius 
r0  =  ctv  /(1  +  z) ,  whereas  r0  might  be  much  smaller.  In  our 
scenario,  both  the  relativistic  outflow  and  the  shell  on  which 
it  impacts  are  particle  energy  dominated;  only  the  shocked 
fluid  is  highly  magnetized.  The  GBM  emission  could  be 
nonthermal  electron  synchrotron  radiation  that  includes 
photospheric  emission  made  at  smaller  radii  than  the 
nonthermal  synchrotron  radiation.  Compton  cooling  of 
shocked  electrons  is  suppressed  in  highly  magnetized 
shocked  plasma,  but  could  be  important  in  a  leptonic  model 
with  external  cocoon  photons  [36].  Our  model  for  GRB 
080916C  also  applies  to  other  Fermi  LAT  GRBs  with 
comparable  T  factors  and  small  beaming  angles. 

Besides  GRBs,  UHECRs  might  also  be  accelerated  in 
other  systems  of  relativistic  outflows,  including  low 
luminosity  GRBs  [37,  38],  radio  galaxies  [39]  and  blazars. 
UHECRs  could  be  formed  through  neutron  escape  when 
photopion  processes  are  important,  which  will  require 
IceCube  neutrino  detections  [40]  to  establish.  In  GRB 
0809 16C,  where  multi-GeV  radiation  is  observed  with 
Fermi,  we  have  shown  that  synchrotron  radiation  from  ultra- 
high  energy  protons  accelerated  in  GRB  jets  explains  the 
delay  of  the  >100  MeV  LAT  emission  with  respect  to  the 
keV  -  MeV  GBM  emission,  and  that  long  duration  GRBs  are 
possible  sites  of  UHECR  acceleration. 
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